THEORETICAL CONSTRAINTS IN THE DESIGN 
OF MULTIVARIABLE CONTROL SYSTEMS 

Final Report 
NAG-1-1361 


January 1993 


Prepared for: 

NASA Langley Research Center 
Hampton, VA 


Prepared by: 
E.G. Rynaski 
D.J. Mook 


CALSPAN-UB RESEARCH CENTER P.O. BOX 400, BUFFALO, NEW YORK 14225 TEL (716) 631 -6900 


TABLE OF CONTENTS 


tkclhm Essz 

Abstract 1 

1.0 Multivariable Control System Design 2 

1.1 Introduction and Background 2 

1.2 Flying Qualities Requirements 3 

2.0 Poles, Zeros and Decoupling 3 

2.1 Introduction 3 

2.2 Transmission Zeros and Decoupling 4 

2.3 Transmission and Transfer Function Zeros 6 

3.0 Dynamic Inversion, Model Following and Direct Design 9 

3.1 Introduction 9 

3.2 Model Following and Dynamic Inversion 10 

3.2. 1 Implicit Model Following 1 1 

3.2.2 Explicit Model Following 12 

3.2.3 General Solution 15 

3.3 Direct Design Methods 19 

4.0 References 24 

Appendix A 

Flight Simulation Derivation A-l 


LIST OF FIGURES AND TABLES 


Figure £ Engs. 

1 Conceptual 2 Input, 3 Output System 7 

2 Implicit Model Following 12 

3 Explicit Model Following 12 

4 Alternate Explicit Model Following Architecture 13 

5 Single Controller Model Following 15 

6 General Model Following Control System 16 

7 Robust Model Following Architecture 18 

Appendix A 

Figure # Page. 

1 Coordinate Axis System A-2 

2 Body Axis System A- 3 

3 Vector Components A- 5 

4 Rotation from Earth-Fixed to Body-Fixed Frame A-8 

5 Stability Axis A-9 

6 Relationship Between the Earth, Body and Stability Axis Frame .... A- 1 0 

7 Short Period and Phugoid Modes A- 17 


Table! Erne 

1 Flight Dynamic Equations A- 13 


ABSTRACT 


The objective of the research performed under this grant was to define and 
investigate the theoretical constraints inherent in the design of multivariable control 
systems. These constraints are manifested by the system transmission zeros that limit or 
bound the areas in which closed loop poles and individual transfer function zeros may be 
placed. These constraints were investigated primarily in the context of system 
decoupling or non-interaction. 

It was proven that decoupling requires the placement of closed loop poles at the 
system transmission zeros. Therefore, the system transmission zeros must be minimum 
phase to guarantee a stable decoupled system. Once decoupling has been accomplished, 
the remaining part of the system exhibits transmission zeros at infinity, so nearly 
complete design freedom is possible in terms of placing both poles and zeros of 
individual closed loop transfer functions. A general, dynamic inversion model following 
system architecture was developed that encompasses both the implicit and explicit 
configuration. Robustness properties are developed along with other attributes of this 
type of system. Finally, a direct design is developed for the longitudinal-vertical degrees 
of freedom of aircraft motion to show how a direct lift flap can be used to improve the 
pitch-heave maneuvering coordination for enhanced flying qualities. 


1.0 Multivariable Control System Design 

1.1 Introduction and Background 

Future flight control systems will be designed to enhance the ability of a pilot to 
fly with more precision than is now generally possible. Precise flight vector control will 
be required for nearly every piloting task, from highly accurate approach, flare and 
landing to ordinance delivery, terrain following/avoidance and high angle of attack 
maneuvering. Most new high performance aircraft have incorporated multiple and often 
redundant control effectors, including canard surfaces, direct lift flaps, thrust vectoring 
and split rudders (ailerons) as part of the controller complement. Relatively little research 
has been accomplished to define how to best use these new and often novel means of 
generating forces and moments on the airplane to enhance the ability of the pilot to fly 
with maximum precision, ease and confidence. 

With only one or two exceptions, aircraft flight control laws are designed on a 
single input-single output basis, and this approach to conceptual design served the 
designer well when the aircraft possessed a conventional controller complement of 
elevator, rudder and aileron. Because additional control effectors can provide for 
enhanced piloting control capability, equivalent methods of design for multicontroller 
configurations are needed now to provide the industry with the conceptual design tools to 
produce superior, multicontroller flight control system configurations. 

Using conventional control surfaces, some maneuvers of an airplane are 
performed indirectly and are possible only because of the aerodynamics coupling inherent 
in the conventional geometry of an airplane. A typical example involves simple vertical 
plane motions - a pilot must normally rotate or pitch an airplane by elevator deflection 
first before an angle of attack is generated that produces a normal force enabling the 
airplane to maneuver (i.e., change flight path). The dynamic delay between pitching and 
normal force generation of a heavy, severely swept wing airplane can make it difficult for 
a pilot to fly (i.e., manipulate his flight path) with precision. A suitable pitch-heavy 
harmony must be provided. On the other hand, many aerodynamic coupling effects 
degrade the ability of a pilot to fly with precision. Aerodynamics coupling produces 
undesirable effects such as adverse/proverse yaw, dutch roll residues in the rolling 
response of an airplane to a pilot stick command and uncoordinated turns. The 
elimination of these undesirable effects often leads to decoupling or non-interaction 
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among the vehicle degrees of freedom of motion. Techniques of multivariable control to 
enhance coupling when desired and decouple when advantageous should be developed to 
a much greater degree than now available. 

1.2 Flying Qualities Requirements 

Flying qualities requirements define the satisfactory and acceptable range of static 
and dynamic response of an airplane to either a pilot command or to environmental 
disturbances. These dynamics are historically expressed in terms of frequency domain 
parameter ranges, such as poles and zeros of particular transfer functions with respect to 
stick and throttle command inputs by the pilot. 

In general, all the poles of a coupled multiple degree of freedom system such as 
an airplane can be exactly specified, or "placed", using state feedback to a single 
controller. But none of the transfer function zeros can be altered, which generally means 
that the aerodynamics coupling among degrees of freedom of motion on the vehicle are 
not significantly changed by the feedback. If the airplane possesses as many independent 
control effectors as degrees of freedom of motion, then all poles and zeros can be exactly 
specified and any closed loop dynamic behavior is possible. The aerodynamic coupling 
among degrees of freedom of motion can be eliminated or altered at will. If the airplane 
possesses more than one control effector but fewer effectors than degrees of freedom of 
motion, then all poles and some individual transfer function zeros may be precisely 
specified, and are constrained by the system transmission zeros. The role of the 
transmission zeros is investigated in this research effort, and it is felt that much progress 
has been made toward the objective of better understanding and insight into the problem. 
This understanding is required before flight control system designers will be able to make 
the best use of the multiple control effectors presently being incorporated into new high 
performance aircraft such as the B-2 bomber, NASP (National Aerospace Plane) and the 
HSCT (High Speed Civil Transport). 

2.0 Poles, Zeros and Decoupling 

2.1 Introduction 

A major objective of a flight control system design is to produce the kind of 
dynamic behavior that will enable the pilot to fly with precision and confidence. In order 
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to accomplish this design, some of the dynamics that produce undesirable aerodynamic 
coupling among degrees of freedom must be removed or the affects eliminated while at 
the same time other dynamic coupling effects should be enhanced or altered. These 
objectives should be attained without significantly increasing the order of the response of 
the aircraft through the introduction of excessive filters and compensation networks and 
without introducing excessive time delay. 

Flying qualities requirements are most often stated in terms of frequency domain 
parameters, such as poles and zeros. This implies that the dynamic behavior of the 
vehicle, to relatively small perturbations from trim, can be described as a linear system. 
Because a linearized representation of the aircraft accurately represents the vehicle 
dynamics, the flying qualities requirements are cast mainly in linear system terms and 
flight control systems are successfully designed on this basis. The analysis in this report 
is applicable to linearized representations of the airplane, although some of the theory, 
such as dynamic inversion, is not restricted to a linear system representation. 

2.2 Transmission Zeros and Decoupling 

The zeros of a single-input, single-output transfer function constrain the area or 
domain of pole placement because the root locus plot involving feedback from the output 
is well ordered. Open loop poles will terminate at the zero locations in the s plane. 

Given the linearized vehicle representation: 

x(t) = Ax(t) + Bu(t) (1) 

or (Is-A) x(s) = Bu(s) in the Laplace domain (2) 

then the numerator polynomial of a transfer function is obtained using Camers' rule. A 
column of the control effectiveness matrix B is substituted for the appropriate column of 


the matrix (Is-A). 

The 

determinant of the resulting matrix yields the 

numerator 

polynomial, i.e.. 






bn 

an... 

ain 




b2i 

S-a22 


= 0 

(3) 
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In the above example, the substitution of bn for the term s-an means that the 
numerator of the Nn(s) transfer function will, in general be represented by a polynomial 
whose order is one less than the order of the denominator defined by D(s) = Ils-Al. 

Cramer's rule can be expanded to obtain transmission zeros for multivariable as 
well as single input, single output systems. Partition equation 4 as shown below in such a 
way that B 2 is non-singular, i.e. IB 2 I 0. 


Xl(t) 


' All 

A12' 

’ Xl(t) ' 

-f 

B, 

. x 2 (t) . 


. A21 

A22 . 

. X2(t) . 


lb 2 J 


u(t) 


(4) 


Cramer's rule can then be used to obtain the transmission zeros, which are the 
singularities of the polynomial 


Is- An Bi 

A21 B2 


(5) 


In general, the order of this polynomial is reduced from D(s) by the number of 
independent inputs to the system, i.e., if B is an n X 2 matrix, the order of the 
transmission zero polynomial will in general be n-2. If there are as many independent 
controllers as degrees of freedom of motion represented by Equation 1, then the system 
will have no finite transmission zeros. 

Using Gauss' algorithm, the determinant of Equation 5 can be reduced to the 
expression given below: 


Is- An Bi 

A21 B2 


^Ils-An+BiB^iM 


( 6 ) 


A control law that will decouple xi(t) from X 2 (t) in Equation 1 is given by 
u(t) = -B2' 1 A 21 Xl(t) + u c (t). 
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The application of this control law to the system of Equation 1 yields 


Xl(t) 


(Ah-BiB2 1 A 2 i) A 12 

Xl ^ 1 

Bi' 

X2(t) . 


0 A 22 

. *2(0 J 

1 B 2 J 


u(t) 


( 7 ) 


which shows that xi(t) is decoupled from X2(t), i.e., motion of xi(t) will have no effect on 
the response X2(t). The eigenvalues of the system of Equation 7 are given by the roots of 
the polynomial 


Ils-A 1 1 +Bi B2' 1 A2ll IIs-A22l = 0 


( 8 ) 


A comparison of Equation 6 and Equation 8 shows that a control law that 
decouples xi(t) from X 2 (t) places poles at the transmission zero locations of the system. 
Because these transmission zeros are invariant with state feedback, they, along with lls- 
A 22 I = 0 define the stability of a closed loop decoupled system. In order to produce a 
decoupled flight control system design, the transmission zeros must be minimum phase. 

It appears that unstable decoupled behavior will not generally be experienced in 
an airplane with conventional geometry and control effectors as long as the control 
surfaces used do not contradict their primary purpose. For instance, to decouple pitch and 
heave degrees of freedom of motion from speed changes using devices designed to 
produce pitching moments and direct lift (i.e., elevator and direct lift flap) would not 
likely produce an unstable speed change dynamics. However, if 0(t) and a(t) decoupling 
were attempted using an elevator and throttle, non-minimum phase transmission zeros are 
possible and perhaps likely. 

2.3 Transmission and Tran sfer Function Zeros 

The transfer function zeros of a single input system are invariant, unaffected by 
feedback. In a multiple input, multiple output system the transmission zeros directly 
affect the zeros of individual closed loop transfer functions. Because the flying qualities 
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of an airplane are strongly influenced by the zeros of individual transfer functions, design 
tools for both pole and zero placement should be developed. 


Consider the two input, three output system depicted schematically in Figure 1 

below: 





Figure 1 Conceptual 2 Input, 3 Output System 


where 


and 


yj = ith output, or quantity to be fed back 
uj= the jth input 


Hi 


the open loop transfer function relating the i^ 1 output to the j^ 1 input 


kji = the feedback function to the j^ 1 input from the i 1 * 1 output 

Nils) 


^J 1 
tij = 


D(s) 


where N(s) are the numerators of the open loop 


transferfunctions and D(s) is the open loop characteristic polynomial 


Using the rules of signal flow developed by Mason, the closed loop transfer 
function relating any output to any input can be written by inspection. For this system 


! 


iE3 


7 


X 


they are 

yi , x_ Nn(s)+k22Li(s)+k23L2(s) 

^ A(s) 

yi , ■ )= Ni2(s)-ki 2 Li(s)-ki3L 2 (s) 
u 2 A(s) 

y2 / 3 -._ N2i(s)+k2iLi(s)+k23L3(s) 

U 1 S " A(s) 

y2 /„\_ ^22( s ) + ki iLi(s)-ki3L3(s) 

U 2 ' A(s) 

y3 ( 3)= N 3 i(s)-k 22 L 3 (s)-k 2 iL 2 (s) 

Ul A(s) 

y 3 , x_ N32(s)+ki iL2(s)+ki3L3(s) 

u 2 ‘ A(s) 


where A(s), the closed loop characteristic polynomial is given by: 


A(s)=D(s)+I kjj Nj i (s)+(k 11 k 2 2 -ki 2 k 2 i)L 1 (s) 

+ Ck 1 1 k 23 -k 1 3 k 2 1 )L 2 (s)+(k 2 2 k 1 3 _ ki2k2 3 )L3(s) ( ^ 0) 

and Li(s), L 2 (s) and L 3 (s) are the system transmission zeros, obtained from 


Li(s)= 


Nn(s)N 22 (s)-N 1 2 (s)N 2 i(s) _D(s)Li(s) 


L 2 (s)= 


D(s) 

N 1 i(s)N 3 2(s)-N 1 2(s)N 31 (s) 


D(s) 


L 3 (s)= 


D(s) 

N 2 1(S)N32(S)-N 13 (S)N22(S) 


(ID 


D(s) 


The minors Nn(s) N 22 (s) - Nj 2 (s) N 2 i(s) must contain the characteristic 
polynomial as a factor (ref. 1). 

From the expressions given above, it can be seen that the numerators of each of 
the transfer functions are directly affected by two of the three transmission zero 
polynomials of the system, so some closed loop transfer function zero placement is 
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possible. Because the order of a transmission zero polynomial such as Li(s) is in general 
one less than an open loop numerator polynomial, the closed loop numerator polynomial 
can be completely specified if n-3 separate transmission zero polynomials are generated 
by feedback. 

It is clear that the closed loop transfer function zeros, such as defined by Ni i(s) + 
k22 Li(s) + k23 L 3 (s) can be obtained from root locus plots involving k22 and k23, as 


k22Ll(s) 

N„(s) 


and 0=1+ 


k23L3(s) 

Nn(s) 


( 12 ) 


where Li(s) and L3(s) represent the invariant root locus plot zeros and Nn(s) represent 
the singularities that will change as a function of the gains k22 and k 23 - A direct design 
approach, in which the closed loop poles and zeros are selected apriori is subject to the 
following constraints 


Nii - N 2 2 - -N 1 2 <1 N2i < r^(s)L 1 (s) 

N 1 1 ciN 3 2 c f N i 2 c 1 N 3 i cl =A(s) L 2 (s) (13) 

N 2 i £l N 3 2 < rN 3 , cl N 22 CI =A(s) L 3 (s) 


where A(s) is the design-objective closed loop characteristic polynomial. 

Because the design method descibed above is complex, and because as many 
design criteria or objectives involve decoupling, the emphasis on design development 
methods is on decoupling, which is in itself a design objective. Direct design methods 
will be stressed, assuming all transmission zeros of the decoupled part of the system will 
be infinite. If this is so, direct dynamic inversion methods will yield the required control 
laws. 

3.0 Dynamic Inversion, Model Following And Direct Design 

3.1 Introduction 

As described in Section 2.2 above, the application of the control law u(t) = -B2' 1 
A 21 xfit) + uc(t) to the partitioned system 
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( 14 ) 


Xl(t) 


' A u 

A12 

’ Xi(t) 

+ 

Bi' 

. X2(t) . 


. A 2 i 

A22 . 

. X2(t) . 


[b 2 J 


will yield a closed loop system in which X2(t) is decoupled from the rest of the system. 
Applying the decoupling control law yields 


xi(t) 


. X2(t) . 



An-BjB^Asi 

0 


A12 

A22 


’ Xl(t) 

. X2(t) . 



U c (t) 


(15) 


and the decoupled part of the system x 2 = A22 x2(t) + B2 u c (t) can be designed 
completely independently as representing the dynamics of interest that will yield 
optimum level 1 flying qualities. The partitioning was chosen such that IB2I * 0, so direct 
dynamic inversion and "exact" model following design methods are feasible. 

3.2 Model Following and Dynamic Inversion 

Assume that a system can be conceptually defined that yields optimum flying 
qualities, the model is 

y(t) = My(t) + B m u(t) (16) 


The objective is to force the plant, defined by x 2 = A22 X2(t) + B2 u c (t) to respond 
"exactly" as the model, defined by Equation 16. The solution is most easily and directly 
obtained by defining a regulator in the error between y(t) and X2(t) as 


e-(A 22 +P) e(t)=0 


(17) 
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where e(t) = y(t)-x 2 (t) (18) 

e(t) = y(t)-x 2 (t) 

and P in any stability matrix selected by the designer. 

Substituting in Equation 17 for e(t) and e(t) yields 

(y(t)-x 2 (t)) - (A 22 +PXy(t)-x 2 (t))=0 (19) 


Substituting one more time for x 2 (t)=A 22 x 2 (t)+B 2 u c (t) yields 


[y(t)-A 22 x 2 (t)-B 2 u c (t)-(A 22 +P)(y(t)-x 2 (t))]=0 (20) 

Finally, solving Equation 20 above for the control input that will force X 2 (t) to 
behave "exactly" as y(t) yields the general control law 

u c (t)=B 2 [y(t)-A 22 y(t)-P(y(t)-x 2 (t))] (21) 

3.2.1 Implicit Model Following 

The general control law of Equation 21 above represents an entire family of 
solutions, depending upon the value of the matrix P selected by the designer, that will 
yield the same dynamic behavior of the plant. If, for instance the matrix P were chosen P 
= M-A 22 , the control law becomes 


u c (t) =B 2 [( M- A 22 )x 2 (t)+B m u m (t)] 


( 22 ) 


The control law indicated by Equation 22 is shown in the block diagram of 
Figure 2 below. The control law involves only feedback of the states x 2 (t). The 
configuration shown in Figure 2 is often referred to as implicit model following. 
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Feedback is used to alter the stability and control derivatives of the plant such that the 
closed loop system behaves as y(t) = My(t)B m u m (t) 



X 2 (t) = y(t) 


Figure 2 Implicit Model Following 
3.2.2 Explicit Model Following 

Since it is well known that almost anything that can be done using feedback can 
also be done in a feedforward sense, an alternate control law canbe obtained that does not 
involve feedback, simply by choosing P=0 in Equation 21. The resulting control law is 
given by 


u c (t)=B 2 [y(t)-A 2 2 y(t)] (23) 

which represents a complete dynamic inversion of the plant, and it can then follow a 
model explicitly. Equation 23 above indicates that y(t) and y(t) are generated 
independently in a computer and the plant controllers are driven through the "gains" B 2 
and - 62^22 by the computer outputs y(t) and y(t). 

The block diagram is shown in Figure 3 below. This architecture is often called 
"explicit" model following 



Figure 3 "Explicit" Model Following 
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As indicated by Equation 23 and shown in Figure 3, the configuration is entirely 
an open loop, feedforward architecture yet yields "exactly" the same results as shown in 
Figure 2, a feedback solution. The control law "gains", B ^ 1 and - 82^22 are a function 
only of the plant, not the model dynamics, so these gain matrices represent a complete 
dynamic inversion of the plant. Because this is so, the model in the computer can be 
changed at will without changing the control law. If, for instance, the flying qualities 
requirements were different for each flying task or flight regime, changing only the model 
dynamics will properly change the dynamic response. If the model represents a 
decoupled vehicle, the plant will respond as a non-interacting system. 

The closed loop or feedback solution shown in Figure 2 and the open loop or 
feedforward solution of Figure 3 represent two extremes of a family of solutions that can 
be obtained that produce equivalent dynamic inversion results, therefore, a complete 
model following capability. There is no control law uniqueness. For instance, consider 
the feedforward or open loop solution of Equation 23. 

u c (t)=B 2 , [y(t)-A 2 2 y(t)] (23) 

which drives the plant through the generation of y(t) and y(t). Two other variations of 
this same control are easily obtained. A second control law is obtained by substituting 
for y(0 in Equation 23 above. 


u c (t)=B 2 I [(M-A 2 2 )y(t)+B m u m (t)] (24) 

This control law strongly involves the computer model dynamics, so as the model 
dynamics are changed, the control law also changes. The block diagram is shown in 
Figure 4 below. 



Figure 4 Alternate "Explicit" Model Following Architecture 
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Because M and A22 are matrices of dimensional stability derivatives and if the 
model is considerably more stable than the plant, i.e., if M » A22 then the solution is 
relatively independent of the plant stability derivatives but still strongly a function of the 
plant control effectiveness terms or stability derivatives. 

A third architecture of the explicit model following control law is given by 
u c (t)=B 2 1 [(l-A 22 M- 1 )y(t)+A 22 M- 1 B m u m (t)] (25) 

This control law reduces simply to u(t)=B2B m u m (t) if the matrix of stability 
derivatives of the plant and the model are identical. 

In summary, because there are three vectors y(t), y(t) and u m (t) of the model and 
any two are required for the model following task, control laws are not unique. In fact, 
feedforward and feedback "hybrid" control law architectures can also be easily obtained. 
Consider the following single controller example given below in which it is desired to 
have an aircraft whose (two degree of freedom) q/de(s) transfer function is given by 

-3-(s) = k(s+1/ W (26) 

6 e s 2 + 2 £cos+co 2 


and it is desired to have the vehicle respond in pitch rate as the model 

km(S+l/'Cm) 

§e S 2 +2^ m C0 m S+t0 m 2 


(27) 


This behavior can be obtained by cascading a lag-lead (or lead-lag) filter model as 
a prefilter, then using pitch rate feedback and an observer to alter the characteristic 
polynomial. The system becomes as shown below in Figure 5 . 
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Figure 5 Single Controller "Model Following" 
where ki and k 2 are obtained from 

s 2 +2£ws+co 2 +k2k(s+l/T e p+kik=s 2 +2£ m G) in s+co m 2 (28) 

or, equating powers of s 

2£ci)+k 2 k+k 1 k=2£ m (D m (29) 

co 2 +k 2 k/T e2 +k i k=0) m 2 

It is clear that although the pitch rate of the aircraft is made to "follow the model" 
using a hybrid feedforward-feedback architecture, the procedure involved characteristics 
of both pitch rate and angle of attack. If this procedure were used to "simulate" the 
dynamics of one airplane using another as the plant, the simulation would not be accurate 
because the more important maneuvering dynamics represented by a(t) would be 
incorrect and would be "higher order," generally detrimental to flying qualities. In 
general, it is not possible to adequately emulate the short period dynamics of a different 
aircraft unless both a pitching moment and direct lift effector were used. 

3.2.3 General Solution 

As indicated in Equation (23) and shown in fig. 3, the feedforward model 
following configuration is entirely an open loop architecture, yet yields "exactly" the 
same result as shown in fig. 2, a feedback solution. The control law "gains", B 2' 1 and 
B2‘*A22 are a function only of the plant, and these gains constitute a complete dynamic 
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inversion of the decoupled part of the plant. Because this is so, the model in the 
computer can be changed at will without changing the control law in any way. If, for 
instance, the flying qualities requirements were different for each flying task or flight 
regime, changing only the model dynamics will properly change the dynamic response. 
In fact, a ground-based simulator operates in exactly the way described above. The only 
difference is that the airplane is moving, not standing still, so the velocity and dynamic 
pressure dependent derivatives A 22 are finite and must be taken into account. 

The closed loop or feedback solution of fig. 2 and the open loop or "feedforward" 
solution of fig. 3 represent two extremes of the family of solutions that can be obtained. 
Neither system is particularly robust in the sense that flight condition variations of the 
stability and control derivatives B 2 and A 22 will cause deviations in the desired response 
unless a lot of gain scheduling is used. Even with gain scheduling, the lack of low 
frequency robustness of the system shown in fig. 3, would preclude its use in an actual 
airplane. However, robustness can be easily provided and in fact, the system shown in 
fig. 3, with feedback added, can result in a robust architecture, accounting for the 
variation of the stability and control derivatives of B 2 and A22- In fact, there is reason to 
believe, as shown below, that minimal gain scheduling is possible for many aircraft. 

The family of solutions that will yield a robust system configuration is given by 
the general solution of Equation (21). The block diagram of the system represented by 
Equation (21) (including the decoupling feedback) is shown in fig. 6 below: 


u c (0 
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y= My+BmUc 

Computer 



Figure 6 General "Model Following" Control System 


16 








--■- v * 





The general solution of Equation (21), (i.e., fig. 5), shows that the matrix P acts on 
the error between the desired output y(t) and the actual output X 2 (t). If B 2 and A 22 are 
known "exactly", the error is zero and P has no function. In fact, P defines not only the 
regulation of the error between y(t) and x2(t) but also defines the perturbation response of 
X 2 (t) to external disturbances (gust alleviation). Most importantly, because P can be an 
integration process, low frequency robustness can be provided, so the response of the 
computer output and the aircraft output will not diverge. Because all the flying qualities 
requirements can be resident in the computer model, gust alleviation or structural mode 
control functions can be done using feedback without affecting the maneuvering flying 
qualities requirements. In fact. Equation 21 is simply a formalization of the present 
common practice of providing feed forward gains in a flight control design. 

If the computer is programmed such that the kinematics of the plant are resident 
in the computer but the aerodynamics in the computer are chosen for flying qualities 
purposes, the computer can represent a global "model following" system. If such a 
system were used for a wide flight range vehicle such as HSCT, the HSCT vehicle would 
always be at the same flight condition as the computer generated "model," but the HSCT 
vehicle would respond dynamically as the "level 1" flying qualities model. There is 
reason to believe that the configuration defined by fig. 6 would work well for a wide 
flight range vehicle. 

Because the function P acts on the error between the computer generated model 
response and the actual vehicle response, P can represent a robust compensator as defined 
by such methods as an H<x> system designed to minimize the error as the stability and 
control derivatives A 22 and B 2 vary with flight condition. 

In general, any compensation network can be expanded in partial fraction 
expansion form to represent proportional, integral and derivative components (or 
designed as a PID system). To show how such a network could be designed for an 
architecture of this type, consider a PID system, resulting in an error control law. 


u e (t) = -K ie (t) - K 2 Je(t)dt - K 3 6(t) 


(30) 
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shown in fig. 7 below: 



Figure 7 Robust Model Following Architecture 

This block diagram is drawn to highlight the effects of the PID feedback and the 
effect of this feedback on the response of the system as stability and control derivatives 
B2 and A22 vary in flight. As shown in the figure, the accuracy of the design depends 
upon the accuracy of B2 and A22, the matrices involved in the dynamic inversion. If the 
feedback gains are made sufficiently large, the system can be made insensitive to 
variations in the stability and control derivatives, i.e., 

if K3 »B2'1 T: insensitive to variations in control effectiveness 

and in control surface nonlinearities 

if Ki »-B2' 1 A22 fc insensitive to variations in B2 1 A22, i-e., the 

stability derivatives 

Making the fair assumption that the feedback is to maintain or improve the 
stability of the vehicle, the integral of the error will guarantee X2(t) = y(t) in the long 
term, regulating always about the trajectory of the airplane on the model computer (for 
instance, a hypersonic NASP with ideal level 1 flying qualities) without the need to 
directly measure some air data (such as velocity) on the vehicle itself. Because the 
integral guarantees long term zero error, the pilot command input-response output is 
independent of flight condition and can be made constant or vary with flight condition 
exactly as specified (Fs vs. n^) in the F.Q. specifications. Because in the long term the 
aircraft steady state response will be totally predictable, a display of the pilot command 
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input will tell the pilot what the steady state flight variable changes will be regardless of 
the sluggishness of the vehicle response, resulting in a "predictor" display. By limiting 
commands at the model, functions such as angle of attack and sideslip can be limited, 
again useful for a NASP system. The problem solving potential of this type of system 
architecture seems significant. 

Even more advantages for this architecture can be realized. For instance, because 

the gains B2' 1 A22 represent a "division" process of dimensional stability and control 

derivatives, many of which are dimensionalized in exactly the same way, many of the 

terms of B2' 1 A22 are simply (constant) ratios of non-dimensional stability derivatives, 
such as c ma /c m( j e , which can be accurately obtained in a wind tunnel. These constants 

occur along the diagonal of the matrix B2‘* A22 , so the robustness oriented gains would 
be designed to minimize the effects of off-diagonal terms, (if their effect is significant). 
Because the aircraft can be made to respond to disturbances entirely independently of the 
computer model (which doesn't respond at all unless turbulence is measured and injected 
into the computer), the feedback can be tailored for gust alleviation or structural mode 
control or even to minimize wing root bending moments (to disturbances) without 
affecting flying qualities (resident in the computer). The minimization of maneuver loads 
can be dealt with in the model computer and reflected in the vehicle itself. For instance, 
if the "model" had a direct lift flap, along with the vehicle, distribution of "model" wing 
loads would be reflected in the vehicle itself. 

3.3 Direct Design Methods 

Because the decoupling described in Section 2.2 was done in such a way that the 
partitioned matrix B2 was non-singular, the remaining system X2(t)=A22 X2(t)+B2Uc(t) 
has no transmission zeros and every element of X2(t) can be exactly "placed," i.e, both 
poles and zeros can be specified. This is shown by an example in this section. 

One of the more important objectives of a Level 1 flight control system design is 
to provide for good pitch/heave harmony in the sense that if the pitch rate and angle of 
attack responses are similar, then the difference between rate of change of attitude and 
rate of change of flight path angle is minimized. The airplane will dynamically "curve" a 
path through the sky. Since d=q+Z a a, then Y=q-a=-Z a a and the rate of change of flight 
path angle is directly proportional to change in angle of attack. Therefore, if the pitch 


19 


rate response is very similar to the angle of attack response, the flight path and pitch 
responses are in harmony. 

Pitch rate and angle of attack vehicle responses are dominated by the two degree 
of freedom short period mode. However, flight precision requires not only attitude/angle 
of attack harmony in the short period mode, but the phugoid mode as well. If the phugoid 
dynamics were decoupled from the pitch rate and angle of attack responses, the 
altitude/flight path angle harmony is maintained with precision beyond the short term 
response of the vehicle. 

Speed change or phugoid dynamics will not affect the angle of attack/pitch rate 
response if velocity changes were decoupled from the rest of the system. If speed 
changes were decoupled, as defined by Equation 7, the subsystem X 2 (t)=A 22 X 2 (t)+B 2 u(t) 
can be represented by the remaining two degree of freedom dynamics 


q(t) 


' MqHx ' 

’ q(t) ' 

4- 

M Se Z 8e 

5 e (t) 

. a(t) . 


1 Za . 

. a(t) . 

T 

. ^8/ Z 5z . 

_5z(t)_ 


(31) 


where de and 5 Z represent deflections of an elevator and direct lift flap (or other direct lift 
force effector). With suitable interconnections to produce force and moment decoupling, 
this equation can be simplified 


q(t) ' 


' Mq M a 

’ q(t) 

4. 

r 

O 

i 

Xa)' 

d(t). 


1 z a 

. a(t) . 


1 

N 

o 

i 

.MO. 


(32) 


The objective will be to define a control law 


u(t)=-k 1 x(t)+L5 st ick(t) 


(33) 
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such that the closed loop system will have a specified closed loop short period frequency 
and damping ratio and the numerators of the a/5 s (s) and q/5 s (s) transfer functions will 
have the same or approximately the same numerator zero. If the feedforward gains L are 
designed such that the pilot commands the elevator and direct lift flap together, the closed 
loop system is defined as 


q(t) 

d(t) 


Mq-Mfekn 

1- Z 5z k 21 


Ma-Mggkn 

' q(t) ' 

_i_ 

Mge 

Z a-Mg z k 2 2. 

. a(t) . 

T 

. Z 8z . 


5 s (t) 


(34) 


The problem is then arranged with four available feedback gains ki i, k22, k 1 2 and 
k 2 i to obtain the short period frequency, damping ratio and numerator zeros desired. The 
closed loop transfer function will have the form 


(35a) 

(35b) 


where and zi, Z 2 , to and £ are to be attained using feedback. 
The closed loop system of Equation 34 is of the form 


a _ Ci(s-zi) 

5 S s 2 +2£cos+co 2 

q_ (s) = C 2 (s-z 2 ) 

5 S s 2 +2£cos+co 2 


a(t) 


’ an ai 2 

' a(t) ' 

J bl l 

. q(t) . 


. a 2 i a 22 . 

. q(t) . 

lb 2 J 


(36) 


whose characteristic polynomial and numerator transfer function polynomials are given 
by 

characteristic polynomial: 

A(s)= s 2 +s(-aii-a 22 )+ana 22 -ai 2 a 2 i 
= s 2 +2£cos+co 2 
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numerator polynomial: 


N a (s) = bi(s-a 2 2+^-ai 2 ) = bi(s-zi) 

bi 



(38a) 

(38b) 


If zi=z 2 , then N^s)=-Z a N a (s)=-Z a bi(s-z 1 ), so the dynamic response of 7 and q differs 
only by a constant. 

From Equation 36, each of the closed loop elements an, an, a2i, and 322 can be 
obtained in terms of the specified poles, zeros and control effectiveness terms as 


a 2 2 


^2^,+z.z, _ W22 k22= w 


(39a) 


an = a 2 2 _ 2^0) - M q -Mg c k22 kn 


M q -ai 1 

kii = — 




(39b) 


an = ^-(a22-zi) = Ma-Mfckn 

b 2 




(39c) 


a 2 i =^-(an-Z2) = 1-Z^k 2 i 

bi 


k 21 =^- 


-8 Z 


(39d) 


The resulting configurations is exactly in the form of an implicit model following 
solution. The only constraint is that zi cannot equal z 2 because element a22 would be 
required to be infinitely large. Other than this, there appears to be no reason why 
complete pole-zero placement cannot be achieved and design criteria attained. 

The feedback methods used are state feedback methods that directly alter stability 
derivatives. Therefore, the order of the system open loop and closed loop is theoretically 
the same and no higher order dynamics or time delay has been deliberately added to the 
system, therefore having no detrimental effect on the flying qualities. Because the 


22 


stability and control derivatives can be directly related to the frequency domain format of 
the flying qualities requirements, the technique described above can be used to design 
directly to the flying qualities specifications. 

Direct design techniques to satisfy flying qualities and precision control criteria 
are still in a stage of initial development, although much progress, as shown in this report, 
has been made. The results shown above demonstrate a specific rather than a totally 
general solution. The ultimate goal of the research is to provide for a completely general 
solution that can be extended to vehicles having many degrees of freedom of motion, 
such as aeroelastic well as rigid body vehicle motions. Although the signal flow oriented 
methods of Section 2 could have been successfully used to accomplish the results shown 
above, it is felt that the complexity of the signal flow methods warrants development of 
direct design methods. 
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APPENDIX A 

Flight Simulation Derivation 

The theoretical developments of the transmission zero work performed under this agreement 
were tested on a simulation of pitch-heave maneuvering coordination for enhanced flying 
qualities. This appendix explains the development of the simulation used in the numerical testing. 

The computer simulations represent the experimental Total In Flight Simulator (TTFS) 
aircraft, operated by the CALSPAN Flight Research Department. 

For simulation purposes, the TIFS is assumed to be flying as a rigid body in a vertical 
(longitudinal) plane, so that there are three degrees-of-freedom (horizontal translation x, vertical 
translation z, and pitch angle 6 ). Aerodynamic coefficients are consistent with CALSPAN reports. 
The details of the simulation follow. 

Coordinate Systems 

In developing the aircraft equations of motion, two different coordinate systems are used, 
shown in Figure 1. One is an inertial coordinate system, which is fixed to the earth and is 
considered to be a non-rotating system. This is a valid assumption since the rotation of the 
earth is negligible in most aircraft dynamic problems. The other system is fixed to the aircraft 
center of gravity and rotates along with the aircraft. This system is referred to as the body axis 
coordinate system. 
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Figure 1 Coordinate Axis System 


Equations of Motion 

Newton’s second law is used to derive the rigid body equations of motion, i.e., the 
conservation of both linear and angular momentum (Nelson [2]): 

= M (1) 

V M = -pH = -j- V f x mv (2) 

^ dt dt^ 

The airplane is considered as a continuum of mass particles (5m) and each elemental mass 
has a velocity ( v) relative to an inertial frame. Newton’s second law is: 







From Figure 2, the velocity of the differential mass can be expressed as: 

dr 

t = v- c + - (4) 

Since the mass of the aircraft is assumed constant and the total mass of the aircraft is found 
by summing the mass elements, Newton’s second law may be written as: 

£^ = ™f + |rS>>“ (5) 

The right hand side of Equation (5) equals zero since r is taken from the center of mass. 
Equation (5) now reduces to: 


F = 



( 6 ) 


The moment equation is developed in terms of the relative velocity of a mass element to 
the center of mass: 


v = v c + — = v c + w x r (7) 

at 

where u is the angular velocity of the aircraft. Substituting this expression into the moment 
Equation (2), the total moment of momentum is: 
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y] (r)8m x Vc *4“ [r* x (a; X r)\8m 


( 8 ) 


H = 


Again, the left hand side of Equation (8) equals zero and the moment of momentum equation 
reduces to: 


H = [fx(wx r)\8m (9) 

A difficulty occurs if the reference is rotating, since the moment will vary with time. To 
eliminate this difficulty, a transformation of the reference frame is made to the body axis system 
(onto the aircraft). A vector A is transformed from a fixed system to a rotating coordinate 
system by: 

dA \ , dA\ . r_ 

it / fixei= ii /’•o*- + b H ( 10 ) 

The force and moment equations, transformed into the body axis system, now become: 


F = m-~~- + m(uJ X v c ) 
at 

(ii) 

X 

it 

(12) 


Vector Components and Scalar Equations 

To obtain the time history results of Equations (11) and (12), it is necessary to express the 
vector equations into component form. The component forms of the forces, moments, gravity, 
linear and angular velocities are shown in Figure 3. 1 The corresponding equations are: 

F = F x i + Fyj + F z k 
M = Li + Mj + Nk 

9 = 9xi + 9 y J + 9zk ^ 

u> = Pi + Qj + Rk 

V p = Ui + Vj + Wk 

-* -* -* 
r — xi + yj + zk 

1 Positive sense is in the direction of the arrows 
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Equation (11) can now be expanded into scalar form: 


F x + mg x =m(u -VR + WQ) 

F y -f mg y = m ^ F + UR — WP'j (14) 

F z + mg z =m(w -UQ + VP^j 







\ 


i E ? 


t -3 


i li 

; c a 

! y 
i y 

! 13 

I e 13 


ff, = (y 2 + !’)Sm - Q^2(xy)Sm- R^(xz)6m 

«y = -PT, (xy)6m + Q^(x 2 z 2 )Sm- R^2(yz)6m ( 15 ) 

//* = — P ^ (xz)5m — Q ^ (yz)£m + (z 2 + y 2 )Sm 

These equations can be expressed in terms of mass moments of inertia about the x, y, and 
z axes: 


Hz = IxxP — IxyQ — IxzR 

Hy = -I xy P + IyyQ ~ Iy Z R ( 16 ) 

H z = -I XZ P - I yz Q + I ZZ R 

Since, for most airplanes, the X-Z plane is a plane of symmetry, the moments of inertia for 
I xy = Iy X = zero. With this assumption, and applying Equations ( 16 ) to the moment expression 
in Equation ( 12 ), the scalar moment equations can be written as: 

L = IxxP - Ix Z R - IxxPQ + {Izz - Iyy)RQ 

M = IyyQ + {Izz - Izz)PR + I*z(P 2 - R Z ) ( 17 ) 

N = I ZZ R - IzzP + {Iyy - Ixx)PQ + hzQR 

where L, M and N are the moments about the X, Y and Z axes, respectively. 

Earth Fixed System and the Kinematic Equations 

The force and moment equations are derived in the body-fixed axis system. But, the position 
of the aircraft must be described by the earth-fixed coordinate system. This is accomplished by 
three consecutive rotations (whose order is important). From Figure 4 , the following rotations 
are made: 

1 . Rotate the X\ Y\Z\ body coordinate frame about the Z\ axis over the yaw angle ($) to the 
X2Y2Z2 body coordinate frame. 

2 . Rotate the X2Y2Z2 body coordinate frame about the Y 2 axis over the pitch angle ( 0 ) to the 
X3Y3Z3 body coordinate frame. 

3 . Rotate the X3Y2Z3 body coordinate frame about the X3 axis over the roll angle ($) to the 
XYZ inertial coordinate frame. 


i J 
i j 
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These angles (yaw, pitch and roll) are known as the Euler angles. 

The velocity components between the body-fixed coordinate system and the earth-fixed 
coordinate system are related by a set of orthogonal transformations (a more complete derivation 
may be found in [1] and [ 2 ]). These transformations are: 



C Q C* 

S$SqC\$ — C$S\$ 

C$S<$Cqi + S$S\& 

m 


Cq 5\j> 

S^SqS^ + C$C\& 

C^SqS^ — S$C\g 

v 

(18) 

-5© 

S$Co 

C*C e 

1 w j 



where (7$ = cos($), = sin('J'), etc. 

The kinematic equations relate the Euler angles ('P, 0 and $) and the angular velocities (P, 
Q and R). From Figure 4, u; must equal the vector sum of the time rate of change of the Euler 
angles about the ki, j’ 2 . and 13 axes: 


U = * k 2 + 0J3 + & (19) 

Using transformations similar to Equation (18) and from Equation (13), the angular velocity 
may be written as: 


UJ 


= Pi + Qj + Rk = i $ sin 0 + + 

j ^ cos 0 sin $ + 0 cos 




0 


k ^ cos 0 cos $ — 0 sin $ 
Equating components, the kinematic equations become (in matrix form): 


( 20 ) 


(P) 


' 1 0 — sin 0 1 



Q 

► = i 

0 cos $ cos 0 sin $ 

> < 

0 > 

UJ 


0 — sin $ cos 0 cos $ J 




( 21 ) 


The time histories of the Euler angles are determined by inverting the 3x3 matrix in Equation 


( 21 ): 


= P + Q sin $ tan 0 + R cos $ tan © 
0 = Q cos $ — R sin $ 

$ = (Q sin $ + R cos $) sec 0 
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( 22 ) 



Aerodynamic Nomenclature 

All forces and moments developed are based on the stability axes system defined by Figure 
5. This system is utilized since experimental data from wind tunnel tests are presented in the 
stability axes. The stability axes are obtained by rotating the body axes ( XY Z ) about the Y = Y S 


axis. This is done over a rotational angle a (known as the angle of attack) until the body-fixed 

— # 

axis (X-axis) coincides with the free stream velocity vector ( V pi ). 



The equation for the angle of attack is: 


a = tan 


-l 


W 

u 


( 23 ) 


The equation for the airspeed is: 

V p = yf(U 2 + V 2 + W 2 ) 


( 24 ) 


The flight path angle (7) is the difference between the pitch angle and the angle of attack 
(7 = 0 — a). Figure 6 depicts the relationship between the different coordinate axes in the 
longitudinal plane. 
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Figure 6 Relationship Between the Earth, Body and Stability Axis Frame 


Forces and Moments 

To simplify the analysis, all forces and moments are developed in the stability axis system. 
The force (F*) and moment ( M a ) components are: 


F 3 = FX/s + F%J, + FX x k> = -Di a + Fyf s - Lk, 
M’ = L A i„ + M A j s + N A k„ 


( 25 ) 


where D = Drag, L = Lift and Fy = Sideforce. L A , M A and N A are aerodynamic moments. 

The steady state forces and moments for a straight line flight are assumed to depend only on 
angle of attack (a), sideslip angle (/?), thrust and the control surface deflections of the elevator 
(5e), ailerons (S A ), rudder (Sr), and aerodynamic coefficients. These dimensionless coefficients 
are comprised of derivatives evaluated at constant Mach and Reynolds number (e.g.): 


Cd — Cd 0 + Ci)q q + Cd, b Se (26) 

where: 

• Cd o = total airplane drag coefficient for a = Se = 0. 

• Coa = total airplane drag change with angle of attack for Se = 0. 

• C Ds b = total airplane drag change with elevator angle for a = 0. 
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A similar analysis is made for the remaining aerodynamic coefficients. 

The aerodynamic forces and moments are expressed in terms of the dimensionless coef- 
ficients, flight dynamic pressure, characteristic length (for the moments only) and a reference 
area: 


D = CnqS = (Cd 0 + Co a a + CD, B ^E S jqS 
F y = Cp Y qS = (c Yf) (3 + C YlA S A + C YfR 8RjqS 
L = C L qS = (C Lo + C La a + C Lf J E )qS 

La = CiqSb 

= (<v + c, lA s A + c K s R + c, p ^ + C,„^)«S6 

M a = CMqSc 

~ (CMq + CMa a + CM( b 8e + CMq^j^jqSc 
Nji — CwqSb 

= (cn p P + Cn, a &a + Cn, r 8r + c Xp2ij- + Cnr 2U]) 


The flight dynamic pressure is : 



where p is the air density and U, is the aircraft’s airspeed. 
The aircraft thrust vector is also divided into components: 

Fj> x = T cos (a) 

Fi= 0 

F^ = — T sin (ct) 

L?p = 0 

My = -Td T 

N} = 0 


(27) 


(28) 


(29) 


(30) 


where dp is the moment arm of thrustline. 

For a trimmed airplane, the forces and moments acting on it are in equilibrium. This 
is accomplished when the pitching moment equals zero and when the lift equals the airplane 
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weight. Using Equations (27) and (28): 


0 = 

^9 r 

~S ^ Ltritn 


(Cm 0 + CM a Ot + Cm^Se') 

(CLo + C La a + C LtB 6 E ) 


(31) 


Solving these equations for a and Se determines the trim value for the angle of attack and 
the trim elevator setting: 


&trim — 


CL trim C MtE - Cl 0 C M(b + Cl, b Cmo 
Cl a C Mi B ~ C Lfj} C Ma 


CL trim C Ma - C La C Ma + C l<i Cm 0 
Cli b Cm„ - C La C MfB 


(32) 

(33) 


Combining the rigid body equations and the aircraft kinematic equations yields a full set 
of flight dynamics equations that describe any aircraft flight path or motion. To accomplish 
this, all forces must first be transformed from the stability axis system to the body axis system. 
From Figure 5: 


F £ = — D cos (a) + L sin (a) + T cos (a) 

F b = F* (34) 

V y 

F b = —D sin (a) — L cos (a) + T sin (a) 

The force equations (14) are then solved in terms of the linear velocity rates (U, V and 
W). Also, the moment equations (17) are solved in terms of the angular velocity rates (P, Q 
and R). To simplify the analysis, let: 

= Ixz{Iyy ~ Ixx) 

h = Izz(hz ~ Iyy) ^ 

h = IxzUzz - Wy) 
h = Ixx(Iyy - Ixx) 

The full flight dynamic motions are now expressed as a 12 th order set of nonlinear differential 
equations (three of which are simple integrations of the linear velocities to yield linear positions). 
These equations can now be integrated to yield flight trajectories, after a transformation is made 
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to the earth-fixed coordinate system by Equation (18), of any variable (e.g. 
airspeed, etc). 

plane altitude, 

U = -gsin (0) + ^ + VR - WQ 

771 

(36) 

F b 

V — g sin ($) cos (0) H — — — UR + WP 

m 

(37) 

W = g cos ($) cos (0) -t — — + UQ — VP 

m 

(38) 

p IxzNa + hz^A - hPQ - I 2 XZ Q R + IxzhzPQ - hRQ 

Ixxhz - l\z 

(39) 

_ M a - (Ixx - Izz)PQ ~ Ixz(P 2 - R 2 ) 
I Y y 

(40) 

p IxxNa + IxzLa + I 2 XZ PQ ~ *3 RQ - h PQ - IxxIxzQR 

Izzlxx - I 2 xz 

(41) 

0 = Q cos $ — R sin $ 

(42) 

^ = (Q sin $ + R cos $) sec 0 

(43) 

$ = P + Q sin $ tan 0 + R cos $ tan 0 

(44) 


Table 1 Flight Dynamic Equations 


The full flight dynamic equations in Table 1 may be linearized into second order differential 
equations by utilizing small disturbance theory. This theory applies to small deviations (for 
angle of attack, sideslip and control surface deflections, etc.) relative to some steady state flight 
condition. For an automatic landing system, this assumption is valid since only small angle 
deviations caused by turbulence or control variables are encountered. This theory is also useful 
in analyzing the stability of an autopilot by using longitudinal transfer functions. 
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Small Perturbation Equations (stability axis) 

Perturbed state equations are derived by replacing all motion variables by a steady state 
value and a perturbation: 

U = U Q + AU V = V 0 + AV 

P = P 0 + AP Q = Q 0 + AQ 

$ = + 0 = 0 O + A© 

This also applies to all forces and moments: 

F x = F Xo + AF X F y = F yo + AFy 

M Mq AM L = Lq A L 


W = W 0 + AW 
R — - Rq -1- AR 
$ = $ 0 -f A$ 


(45) 


F z = F Za + AF Z 
N = No + AN 


(46) 


The change in forces and moments can be expressed in terms of the perturbation variables 
by using a Taylor series expansion: 


A F. = A U + W + A S E + ^A S t 


1 8U 

dF v 


dw 


88e 


d6x 


dFy 

dP 


AFy = ^AV + ^ A P + ^ A R + sJL AS R 


dFy 

dR 


dFy 

df>R 


Ar - - w“ * w" * * §"■ * 


(47) 


A£ = w AV + § AP + i A * + w r as « + w a A6a 


. dN ATr dN 8N dN dN A _ 

= — A^ + W AP + ^AiE + ^A*« + ^Afc 


The partial derivatives in Equation (47) are called the stability derivatives. For convenience, 
these derivatives are divided by the aircraft mass. These new symbols are defined as dimensional 
stability derivatives (e.g.): 


X u = 




(48) 
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Replacing the forces and moments with Equation (25) the longitudinal dimensional stability 



derivatives, in the stability axis system, become: 




y ~2 Cd 0 <1oS 

“ ~ mUo 

-2C DfB q 0 S 

771 

„ -2CL 0 qoS 

u ” mU 0 


; r? 

ry -CL a qoSc 

a 2 mU 0 

-ClqQqSc 

q ~ 2 mU 0 

ry -Cl.JoS 

Zt ‘ = m 


■ — 

,, 2(7^0 9o5 , c 

M “ = lyyUo 

-CM a qoSc 

JVlCL r 

Iyy 

,, Cm^qSc 2 

Ma ~ 2 I Y yUq 

( 49 ) 

; 5 

., CmoQ oSc 2 

M '= 2lyyU„ 

, , CM, B qoSc 

m Sb = 

Iyy 

„ CmtJqSc. 

M T a - J 

Iyy 


! ?! 

x a = d&s. 

- Ci 0 )qoS r , 

Zjql — 
771 

-(Cl,, + CD 0 )qoS 

m 



These dimensional stability derivatives, along with the force and moment perturbation 
equations, can now be evaluated into the full flight dynamic equations (Table 1). To simplify 
the analysis, for small angle deflections, let: 


cos 0 « 1.0 
sin 0 ~ 0 


(50) 


Taking a Laplace Transformation of the resulting equations determines the longitudinal 
transfer functions with the elevator setting (Se) as the input and the horizontal velocity component 
(U), angle of attack (a) and pitch angle (0) as the output variables. These longitudinal transfer 
functions, in matrix form, are: 


r 

A12 

A 13 ' 


, U ( a l \ 

6 b (*) I 

< A21 

A 22 

A 23 

> < 

■ a (/I ( 

,^31 

A 32 

>133 , 


M I 

6 b (») ' 


Xs. ) 

Z Sb 

m Sb J 


( 51 ) 
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where: 


-An = (s — X u ) 

An = -X a 
An = 0cos 0 O 
A21 = —Z u 

An — s{Uq — Za ) — Z a 
An = ~(Z q + Uo)s sin 0o 
An = ~{M U + Mr.) 

A 32 = —(sM& + M a ) 

An = s 2 — 3 Mg 


( 52 ) 


Taking the inverse of Equation (51) yields the transfer functions for the longitudinal mode. 
Using Cramer’s Rule, the denominator of the transfer functions is a fourth order polynomial. 
The roots of this polynomial form the characteristic equation and determine the stability of an 
aircraft. The characteristic equation is represented by two oscillatory modes of motion (short 
period and phugoid). 


Short and Long Period Approximations 

The longitudinal motion of an aircraft is determined by two modes. The first is the short 
period mode. This is characterized by a highly damped, high natural frequency oscillation at 
approximately constant speed (A U « 0). The other mode is the long period or phugoid mode. 
This is caused by a gradual change between potential and kinetic energy and is characterized 
by a low damped, low natural frequency oscillation at approximately constant angle of attack 
(Aa ~ 0). Both modes are determined by the characteristic roots of Equation (51). The concept 
of these modes is illustrated in Figure 7. 
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Figure 7 Short Period and Phugoid Modes 

The short period mode is obtained by approximately neglecting the velocity term in Equation 
(51). This equation now reduces to: 



Applying the inverse to this equation: 

q(s) sZs B + (M 6b U 0 - M q Z 6s ) 

S E {s) ~ s{s 2 U 0 - s{M q U 0 + Z a + U 0 M a ) + ( Z a M q - M a U 0 )} 

( 54 ) 

Q(s) s{UqMs b + Zs B Mg) + (MgZsg — Z a Mj B ) 

S E (») ~ *U 2 Uo - s(M q U 0 + Z a + U 0 Ma) + ( Z a M q - M a U 0 )} 

Comparing the denominator of this equation to the standard frequency form of (s 2 + 
s2(u> nsp + the natural frequency and damping ratio of the short period mode are: 
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! ii nr pk 


( 55 ) 


W »SP 



M a 


, -(M,U„ + Z„ + MaUo) 

Csp - 2^Uo 

The phugoid approximation is derived by neglecting the angle of attack term in Equation 
(51). The resulting equations for this mode become: 


Ujs) _ Z SB Uog 

8e{s) (s*Uo — sX u Uq — Z u g ) 

©(*) = 

S E {s) ( s 2 U 0 - sX u U Q - Z u g ) 


(56) 




The natural frequency and damping ratio for the phugoid mode are: 


v np 



(p - 


-X u 

2cL» np 


(57) 


The combination of the short period and phugoid modes describes the complete longitudinal 
aircraft motion (i.e.): 


QQ>) 

8e( s ) 


kes(T®i± + 1 )(^ r e 2 3 + 1) 




(58) 
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